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Abstract 

We study the validity of a relation by Drell, Levy and Yan (DLY) connecting the deep 
inelastic structure (DIS) functions and the single-particle fragmentation functions in e + e~ 
annihilation which are defined in the spacelike (q 2 < 0) and timelike (q 2 > 0) regions 
respectively. Here q denotes the momentum of the virtual photon exchanged in the deep 
inelastic scattering process or the annihilation process. An extension of the DLY-relation, 
which originally was only derived in the scaling parton model, to all orders in QCD leads 
to a connection between the two evolution kernels determining the g 2 -dependence of the 
DIS structure functions and the fragmentation functions respectively. In relation to this 
we derive the transformation relations between the space-and time-like splitting functions 
up to next-to-leading order (NLO) and the coefficient functions up to NNLO both for 
unpolarized and polarized scattering. It is shown that the evolution kernels describing the 
combined singlet evolution for the structure functions ^(z, Q 2 ), Fl(x, Q 2 ) where Q 2 = \q 2 \ 
or F2(x,Q 2 ),dF2(x,Q 2 )/dln(Q 2 ) and the corresponding fragmentation functions satisfy 
the DLY relation up to next-to- leading order. We also comment on a relation proposed by 
Gribov and Lipatov. 



1 Introduction 



Neutral current deep inelastic lepton-nucleon scattering and single nucleon inclusive production 
in e + e~ pair-annihilation are formally related by crossing the kinematic channels. Already 
before the advent of Quantum Chromodynamics (QCD) Drell, Levy, and Yan mentioned the 
possibility [[[], that the deep inelastic scattering structure functions at the one side and the 
nucleon fragmentation functions in e + e~ pair-annihilation on the other side may be related by 
an analytic continuation from the t- to the s-channel. The hadronic tensors for the space-like 
process of deep inelastic scattering (DIS) and the time-like single nucleon inclusive reaction 
would therefore be related by 

W${q,p) = -WjjP{q,-p) . (1) 

Here p denotes the nucleon momentum and q is the 4-momentum transfer to the hadronic system, 
with q 2 < for deep inelastic scattering and q 2 > for e + e~-annihilation. 

At that time the physical nucleons were considered as bound states built up of bare nucleons 
and pions in the context of the Yukawa theory. The interactions were described by Bethe- 
Salpeter || or Faddeev-type || equations and their generalization || [] aiming at a perturbative 
description of the structure and fragmentation functions. In these theories neither infra-red nor 
collinear singularities are occurring. One may think off a general representation of the structure 
and fragmentation functions in terms of current-current expectation values. However it was 
already shown [p]] that for e + e~-annihilation also diagrams of distinct connectedness appear (cf. 
also |§) which are absent in DIS so that a proof of Eq. (fj) at the non-perturbative level becomes 
very difficult. The relation could be established for the aforementioned ladder-models |I|, [8|, [J Q. 
Within QCD the picture changes. Here it turns out that a thorough perturbative description 
of the structure functions and fragmentation functions is not possible. However, perturbation 
theory may be used to describe the scaling violations of these functions at large values of \q 2 \ 
where the running coupling constant is small. The QCD-improved parton model, moreover, 
exhibits a similarity with the approach by Drell, Levy, and Yan, since in the range where single 
parton states are dominating, i.e. for the contributions of lowest twist, the non-perturbative 
contributions factorize. One therefore may calculate the respective one-particle evolution kernels 
and study their behavior under the crossing from the t- to the s-channel. A further complication 
in the case of a vector-theory as QCD is the emergence of infrared and also collinear singularities 
which have an essential impact on the crossing because of the behavior of the kernels at x — 1 0. 
Here x denotes the Bjorken scaling variable which will be defined differently for the timelike 
and spacelike region except for x = 1 where both definitions lead to the same value for x. As 



a consequence of the Bloch-Nordsiek theorem |L3[ the kernels become distribution-valued for 



x = 1 Ell- 



in leading order for unpolarized scattering the crossing relations mentioned above were given 
in [[Tj|. Similar relations hold in the polarized case. In this order these kernels are nothing but the 
lowest order splitting functions P$\x), which are obtained from the inverse Mellin transforms 



kl 
jV(0) 



of the anomalous dimensions 7&> fLE 



1 Later on massive vector meson ladder-models were studied in [|| Q], where the crossing relation Eq. ([l]) was 
verified for the respective kernels. 

2 For a review on the early developments see [ [To| . 

3 Eq. (|l|) was originally postulated assuming that those terms are absent [Q. See also the subsequent discussion 
in [p"l| pointing out that the exponent of the structure functions oc (1 — x) p near x = 1 H] needs not to be integer. 
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It is the aim of the present paper to investigate the validity of the Drell-Levy-Yan (DLY) 
relation, if applied to perturbatively calculable partonic structure functions and quantities related 
to them up to the level of two-loop order. To establish this crossing relation between space- 
and time-like processes one has to study scheme-invariant quantities which are the physical 
evolution kernels for specific choices of observables as the unpolarized and polarized structure 
and fragmentation functions or derivatives of them w.r.t. q 2 . Furthermore, conditions are derived 
for the transformation of the splitting and coefficient functions from the space- to the time-like 
case. For the coefficient functions we extend the discussion to the NNLO level. Other relations 
between the splitting functions such as supersymmetric relations and relations due to conformal 



symmetry were discussed elsewhere, cf. e.g. [17 



The paper is organized as follows. Basic relations for the deep inelastic structure and frag- 
mentation functions are summarized in section 2. In section 3 scheme-invariant combinations 
of coefficient and splitting functions are constructed for the space- and time-like processes both 
for unpolarized and polarized deep inelastic reactions where we consider two principal examples. 
The Drell-Levy-Yan relation is studied in detail in section 4. We also comment on a relation by 



Gribov and Lipatov |19J which emerged in the same context. Section 5 contains the conclusions. 
In the appendix we present the differences between to the space- and time-like coefficient func- 
tions at 0(a 2 ) as well as the convolution relations which are needed for the investigation of the 
DLY-relation. 

2 Structure Functions and Fragmentation Functions 

Deep inelastic scattering (DIS) of a lepton (I) off a hadron target (P) is described by the process 

1(h) + P(jp) -> l(k 2 ) + 'X' , q = k 1 -k 2 , q 2 = -Q 2 < . (2) 

where l X' represents an inclusive final state. When a single gauge boson is exchanged between 
the incoming lepton and the hadron the above process factorizes into the leptonic part and the 
remaining hadronic part. In the case of forward scattering the scattering matrix element can be 
written in terms of the leptonic tensor L^ u and the hadronic tensor W^ v by 

\M\ 2 = L^W, U . (3) 



The hadronic tensor |20| contains the unpolarized and polarized deep inelastic structure functions 
Fi and g>j. If the process is mediated by photon only we have i = 1, 2 in both the polarized and 
unpolarized case. Notice that instead of F± we can also take the longitudinal structure function 
Fl. At asymptotic values of the kinematic variables structure functions only depend on Q 2 and 
the Bjorken scaling variable 

x B = 7T-, < x B < 1 . (4) 

In QCD the Q 2 dependence of the structure functions is only logarithmic and it accounts for the 
violation of scaling. In the context of the parton model the structure functions can be expressed 
in terms c 
(k = q,g) 



(s) 

in terms of quark and gluon densities and the corresponding spacelike coefficient functions (7 i fe 



4 Similar relations hold for the polarized structure functions gi(x, Q 2 ) and gi(x, Q 2 ) on the level of twist 2. 
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x C, 



( S ) ( z 9l\ + fNs(^ 2 



r (S)NS ( 

*> 9 P/4 



2,L. 



(5) 



where ej denotes the charge of the jth quark flavor and Nf represents the number of light flavors. 
The scale appearing in the above equation, denotes the factorization scale which is introduced 
while removing the collinear singularities from the partonic structure functions. In addition one 
encounters a dependence on the renormalization scale fi r which arises in the renormalization 
procedure. For convenience this scale is put equal to the factorization scale in the following. 
Notice that the structure functions Fj and gi do not depend on these scales. However the parton 
densities and the coefficient functions, which do depend on these scales, satisfy renormalization 
group equations which will be shown below. 

In Eq. (|5|) the index (S) in the structure functions indicates the space-like nature of the 
process (q 2 < 0). Furthermore in Eq. (^) appear the singlet 'S' and non-singlet 'NS' combinations 
of parton densities which are defined by 



E 



u fa /4) + /» (*> $ 



and 



4f M) = f4*>fi) + f*(*>fi)-wf!(*>ti) 



(6) 



(7) 



respectively. Corresponding formulae hold for polarized scattering. In this case the polarized 
parton densities and polarized coefficient functions are denoted by Af k (z,[£f) (k = q,g) and 
AQ, k (z,Q 2 /fi}) (z = l,2). 

Whereas in deep inelastic scattering the constituent structure of the nucleons is studied, 
hadroproduction at e + e~ colliders provides us with information about the fragmentation process 
of these constituents into the hadrons. This information is contained in the fragmentation 
functions observed in the reaction []J 



1(h) + l(k 2 ) P{p) + l X' , q = k 1 + k 2 



q 2 = Q 2 >0 



where the symbols have he same meaning as in Eq. ([|). These fragmentation functions are 
the analogues of the DIS structure functions Therefore in the QCD improved parton model 
these functions can be expressed in a similar way in terms of parton fragmentation densities 
(k = q, g) multiplied by timelike coefficient functions i.e. 



rif j 

FP(x El Q 2 ) = x E J2eU 

i=1 Jx E 



1 dz 



■Di 



Xe 2 



C 



(T)S 




+ D 



x E 

z 



x [z, Qfj + D^ s /ij) CZ WS (z, ^ 



i = 2,L, 



where the corresponding scaling variable for the process in Eq. (9) is defined by 

2p.q 



x E 



CP 



< x E < 1 



(9) 



(10) 
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The symbol T appearing within parentheses in Eq. @ denotes that the fragmentation functions 
are measured in time-like processes. The scales fif and /i r are defined in the same way as in 
Eq. (^) where like in DIS we set the renormalization scale equal to the factorization scale. 
Furthermore the definitions for the singlet and non-singlet parton fragmentation functions are 
the same as those for the parton densities given in Eqs. (|| [7]). Similarly as in DIS one can also 
study the annihilation processes in Eq. (|]) where the hadron P is polarized. This entails the 

(T) (T) 

definition of the polarized fragmentation functions denoted by g\ and g\ for which one can 
present a similar formula as in Eq. (|9|). Very often one also encounters the transverse structure 
function which in the timelike and spacelike case is given by 



1 

2x 



x,Q 2 



F ( L R \x,Q 2 ) 



with 



R = S (x = xb) 



R = T (x 



X E 



'ID 



3 Scheme— invariant Combinations 

In this section we give a short outline of the origin of the factorization scheme dependence 
of the anomalous dimensions (splitting functions) and the coefficient functions. We also show 
how this dependence disappears in the evolution of the structure functions w.r.t. the kinematic 
variable Q 2 . The discussion below deals with the DIS structure functions but the conclusions 
also hold for the fragmentation functions. The partonic structure functions denoted by J-^ 
(i = 1,2, L, k = q,g), representing the QCD radiative corrections, contain various divergences. 
First these divergences have to be regularized for which the most convenient way is to choose the 
method of n-dimensional regularization. Using this method the singularities reveal themselves 
in the form of pole terms of the type 1/e 7 , with n = 4 + e, in the quantity Tij-. The infrared 
divergences cancel between virtual and bremsstrahlung contributions by virtue of the Bloch- 
Nordsieck theorem [13]. Due to the Kinoshita-Lee-Nauenberg theorem |H] all the final state 



mass singularities are canceled too since the DIS structure function is an inclusive quantity. 
Then one is left with only two types of singularities. The first one originates from the ultraviolet 
region. This type of singularities is removed via a redefinition of the parameters appearing in the 
QCD Lagrangian. An example is the coupling constant which becomes equal to a s (n 2 ) where 
fi r is the renormalization scale. After coupling constant renormalization the hadronic structure 
function can be written as follows 

Fi(x, Q 2 ) = £ (f lk (a s (n 2 ), %, ^ e) ® f k ) (x) , (12) 

k=q,g \ V / 

where the symbol <S> denotes the Mellin-convolution defined by 

{f®g){z)=f dzi [ dz 2 f(z 1 )g(z 2 )S(z- z x z 2 ) . (13) 



JO 

Furthermore f k is defined as the bare parton density which is scale independent and is an 
unphysical object because of the singular behavior of T^. Notice that the latter depends on 
the scale /i r and therefore on the renormalization scheme w.r.t. the coupling constant. The 
parameter \x originates from n-dimensional regularization because in this method the coupling 
constant gets a dimension. The second type of singularity originates from the collinear region 
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which can be attributed to the vanishing mass of the initial state parton represented by either 
the (anti-) quark or the gluon. Hence the e in Eq. ( T% represents the collinear singularities 
which are removed from the partonic structure function via mass factorization and transferred 
to a transition function as follows 




3.f|)w.(«.«)|4«)) 



This procedure provides us with the coefficient function denoted by CV 
(0) into Eq. {[lj) leads to the result 



(14) 



--1,9 




Substitution of Eq. 



(15) 



where the renormalized parton density is defined as 



ft(z,a s (i4 



/4 



£ (r lk (a s (tf),^,^,e)®f k )(z) 



2 / 



=q,9 



(16) 



Since the mass factorization can be carried out in various ways one is left with an additional 
scheme dependence which comes on top of the renormalization scheme dependence entering 
the coupling constant in Eq. ([12]). The former only shows up in the parton densities and 
the coefficient functions and it only disappears in specific combinations representing physical 
quantities. Hence physical quantities are invariant under scheme transformation. Like in the 
case of renormalization, mass factorization leads to the introduction of a scale fif called mass 
factorization scale which is related to the factorization scheme dependence. Like in the latter case 
/!/ drops out in physical quantities as the DIS structure functions or fragmentation functions. 
The change of the parton densities and the coefficient functions with respect to a variation 
in the scales fi r and /i/ is determined by the renormalization group equation (RGE) |22| The 
renormalization group equation of the parton densities follow from the one presented for the 
transition functions Ti k . The latter takes the following form 



'9a s (/iV 



V mh (a s {iij), jt, l,e)J (z) 







47T 



6(1 -z) 



(17) 



where we have set /i r — /if for simplicity. The functions Pij{a s ,e, z^j appearing in the above 
equation are the splitting functions. Furthermore the beta-function is defined by 



d $ 



-f3 a 2 s (fi 2 r ) - (3 1 a 3 s (/j 2 '' 



The same equation as in Eq. (|T^) also applies to the parton density because of the definition in 
Eq. flTB"!) . The scale dependence of the coefficient function in Eq. ( Tl5|) is given by 



(9 1 



Cim( a s(^ 2 f)i — 5" ) 1 




0. 



(19) 
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As has been mentioned above scheme transformations such as 

Tik — > Zi m £g> r m fc , Cij — > ^2 A,m ® -^m/ 1 5 (20) 

m=q,g m=q,g 

will not alter the physical observable like e.g. the structure functions and fragmentation func- 
tions. The relation between the splitting and coefficient functions computed in two different 
schemes is found to be 

Pik = Z lm (g) P mn (g) (Z~ 1 ) nk -2f3(a s ) Y z im® -^—(Z~ l ) mk} (21) 

{m,n}=q,g rn=q,g « fl s 

Ci,l = ^2 Ci,m <8> (Z ) m j . (22) 
m=q,g 

Below we present the relation between the coefficient functions computed in two different schemes 
up to order a 2 (Q 2 ). Notice that we have chosen here fi 2 = Q 2 in order to get rid off the logarithms 
\n(Q 2 //ij) which usually appear. Up to 0(a 2 ) one obtains 

Ci, q = 6(1 -z) + a s fc« + Z$) + a 2 s (c® + Zg> + (Z%) 2 + Z« ® Z« 



Ci,g 



+CV®ZU+CV®ZM)+... , (23) 

«- (Cfi? + *$) + «2 (Cff + Z$ + Z$ ® (Z$ + Zg>) 
+CS®Z$+CV® Z $)+.... (24) 



In the subsequent part of the paper it is much more convenient to derive the expressions in the 
Mellin transform representation so that one can avoid the convolution symbol <8>. The Mellin 
transform of a function f(z) is given by 



f {N) = / dz z N ~' f(z) (25) 
Jo 

In this way Eq. ( |1"3"D can be written as 

(/ ® 9) N = t dz z N ~ l (f ® g)(z) = f-g N . (26) 
Jo 

Since the structure functions are scheme independent they become renormalization group invari- 
ants. Hence they satisfy the RG equation 



F l N (x,Q 2 )=0 (27) 



The equation above follows from combining Eqs. flT7| , pT9| ) and fll5|) . However the independence of 
the structure function on the scales fif and \i r is not manifest when multiplying parton densities 
with coefficient functions. In particular when the perturbation series of a physical quantity is 
computed up to finite order there is a residual dependence on these unphysical scales (see e.g. 
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23| . Their influence is expected to become smaller when higher order terms in the perturbation 

To avoid the problem of the factorization scheme dependence 



series are taken into account 
of the structure function when the perturbation series is truncated up to fixed order it is better 
to study evolution equations for the structure functions with respect to a physical scale which is 
represented by a kinematic variable like Q 2 . In these type of evolution equations the kernels are 
factorization scheme independent order by order in perturbation theory. However the dependence 
on the choice of renormalization scheme and therefore the dependence on \i r remains so that 
one is able to obtain a better estimate of the theoretical error on a s . For such an equation one 
needs two different structure functions called Fa(x, Q 2 ) and Fb(x, Q 2 ). Examples are A = 2 and 
B = L or Fa and Q 2 d Fa/cI Q 2 . Limiting ourselves to the singlet case, the evolution equation 
for the non-singlet structure functions is even more simple, one can write 



Ql 



Ql 



^I,g \ a s(f j2 f), 



A,B. 



(28) 



Here one can view the Cf t (I = A,B, I = q,g) as matrix elements so that the equation above 
has the form 

r<N r<N \ f tN 

^Aq U Ag \ | J a 



TpN 
B 



(29) 



The coefficient functions satisfy the RG-equation in Eq. (|T^) and the solution is given by the 
T-ordered exponential |25[ 



Cf, m {a s {Q 2 ),l) [T as 



exp 



a s (Q 2 ) 7 ^( X ) 

dx — — 

a 3 ( M 2) ' 2/3(x) 



ml 



where 7^ is the anomalous dimension matrix defined by 



— / dzz 



N-l 



Pik{z) 



(30) 



(31) 



We will now differentiate the coefficient functions w.r.t. Q 2 



Q' 



dC^ k (a s ^ 2 ),Q 2 /^ 2 ) 
dQ 2 

dC? m {a s {Q 2 ) 



P(a s (Q 2 )) 



(3{a s {Q 2 ) 



(C N 



m,J 



da s (Q 2 ) 

-\cf tm {a s {Q% l) 7 lMQ 2 )) (c*)^ (a s (Q 2 



dC? ik (a.(ji 2 f ),Q 2 /fi) 
da s (Q 2 ) 



(a s (Q 2 )A) 



Q 2 



(32) 



One can show that the expression above is invariant under scheme transformations. The latter 
are given by 



{m,n}=q,g 



m=q,g 



£^(n^> ( cN )u= E z ?™{c N y". 



(33) 
(34) 



m=q,g 



m=q,g 
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Since the Q 2 "dependence only resides in the coefficient function the same evolution equation as 
in Eq. (|32p also applies to Ff in Eq. (ESI). For a short-hand notation we introduce the evolution 
variable t 

2 / a,(Q 2 ) \ 

i = -/% ln UM)J' (35) 

so that we obtain 

d ( Ff \ _ \(Kf A K% \ ( Ff 
dt\F%)--\\ Kg A K% ) { F» 

where the physical (scheme invariant) kernel is given by 

A 9Cj (t) /„jv\-i , , f3 a s (Q 2 ) N N ,,,{ riN \- 1 



(36) 



K N 



dt 



(37) 



The kernels Kfj depend both on the anomalous dimensions 7/^ and the coefficient functions 
Cfi(Q 2 ) but the latter two quantities are combined in a factorization scheme independent way. 
This factorization scheme independence of Kfj holds order by order in perturbation theory. 
Using the series expansions for the anomalous dimensions and coefficient functions in terms of 
the strong coupling constant 

00 i , 00 , , 

Tfi = E < + \Q 2 ) (7*9 1 ' C^(Q 2 ) = j:<(Q 2 )(C N )^ , l,k = q,g , I = A,B, 

n=0 ra=0 

(38) 

one can compute order by order the coefficients in the perturbation series of the kernel 

00 , > 

Kfj = Y.<{Q 2 ){K N tj ■ ( 39 ) 



n=0 



Notice that the coefficients [K J ^ are not invariant with respect to a finite renormalization of 
the coupling constant. This dependence is removed when the perturbation series in Eq. ( |3"§D is 
resummed in all orders. 



3.1 F 2 (x, Q 2 ) and F L (x, Q 2 ) 

Let us consider now two specific examples, choosing the structure functions F 2 (x,Q 2 ) and 
Fl(x,Q 2 ) or the structure function F 2 (x,Q 2 ) and its slope dF 2 (x,Q 2 )/dt as the observables 
Fa,b{x, Q 2 )- in this combination of observables it is convenient to normalize the structure func- 
tion F L (x,Q 2 ) to its gluonic contribution in lowest order. This is because F L (x,Q 2 ) vanishes 
in zeroth order of a s due to the Callan-Gross relation, cf. Eq. QTTp, contrary to the structure 
function F 2 (x, Q 2 ). Therefore this normalization accounts for keeping the same order in the 
coupling constant for the two quantities 

f a(Q 2 ) = if (5) (Q 2 ), F-(Q 2 ) = f£ ( 5i (1) • (40) 

Since both the coefficient functions C£^ and Cf^ are scheme invariants such a normalization is 
possible. We expand now the kernels Kfj for this choice of observables into a series in a s . The 
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lowest order contribution is well-known, cf. e.g. 

C N(1) 

c 



K N(0) _ N{0 ) _ ^L,q N(0) R N(0) _ JV(0) 

^22 ~ lqq n N{l)^qg n 2L ~ Iqg 



'L,g 



' r N(l)\ 2 r N(l) 
K N(0) _ N(0) _ t ^L,g N(0) K N(0) _ N{0) , U L,g JV(Q) 

IV L2 Igg „N(\) Iqg J ^ LL Inn ^,Nf~\) Iqg 



i gq 

,N(1) 



n N{l) I iqg LL Igg 1 „JV(1) '99 

U L,g ) U L,g 



I U L,q / N(0) _ n Jf(f))\ 

^^N^ytqq hg J 
u L,g 



To 



next-to-leading order in 

( JV(2) 



a s(Q 2 ), one finds 

a r N(l) ( 

K N M - -MP ft,,w(Q) °£,g ( 
-2 "/•/ "/'/ qN(i) 1 • 



K. 



,iV(l) 

gg 



_ Pi MO) \ , U L,g n N(l) N(0) 

'gg J ^ c N(i) u 2,9 r OT 



^L,q I ^L,q 

r N(l) " r 



°i,9 



^(1) 



/So 

~/V(l) ^(2)' 

°L,q U L,g 

°i,8 °I,9 



N(0) + qN(1) N(0) 
Iqg ^ ^2,g I gq 



r N(l) ( 
^L,q r N(l) V(Q) . r JV(l) 

JV(i) ^2,9 7 99 + ^Po ^2,g 



^iV(l) 

r AT(l) 
~,AT(1) ^2,9 



L,g 



Ml) _ . _. 

21. Iqg Qlqg 1 2,y 

^(1) 



^i^(o) _ ^(l)/ jv(o) _ tv(o)n , o/9nr 7V(1) 

^ Iqg ^2,g {lqq Igg ) + ^Po^2, S 

Po 



^L,q n N(l) ^L,g 
^777^9 n KFT- 



C N(1) 



K N{1) ~ T 
A L2 - 7 OT 



Pi AT(0) , 



+ 



mi) _ ^.M) 



C L,q ] ( Ml) _ &N(0) S 

Igg lgg 



,N(1) 



n N(2) n N(l) 

U L,g _ )^q_ r N{l) --, , 

n N{l) r JV(l)°2>9 n N{l) 

°i,9 V°i,9 / 



r 7V(l)\ 2 

I ^(1) 



ryMO) 

lqq 



n N(l) \ 3 n N(l) n N{2) 



h\ < (1) + 2 



r N(l) 



n N{\) n N{l) 
^L,g ^L,g 



qN(1) \ 2 ^AT(2) 



'L,g 

n N(l) I ^(1) 
/ °i,9 



c 



r N{l) 
°i,9 



^JV(l) 
L/ 2,g 



'qg ^ 



C 



N(l) 
L,q 



c 



a 



N(l) ^2,9 
L,g 



JV(1) 



^(1) , 
°2,q + 



r N{2) \ 
^ L >9 ^JV(O) 
N(l) I '9<7 



i,9 / 



10 



' r N{2) 




r N(l) " t " 





^JV(l) 



U L,q r N(l) 
N{1) °2,g 



c 



L,g 



1 99 



' n N(2) n N{l) n N{2)\ 
^L,q °r " u ' ■ « l 



^AT(l) 



jag ^&,g 

r AT(l) ^(1) 



(44) 



jV(l) 
LL 



N(l) _ Pi N(0) , ^L,g N (l) 
1 99 p Q 1 99 ~ qN(1) y ! Q9 



Pi N(0) 



U L,q 
' r N(l) 



^2,g Iqq + 



r N{2) 
^L,q 

r N(l) 



n N(l) n N(2) 

U L,q U L,g 

n N{l) n N(l) 

U L,g ^L,g 



| U L,g r N(l) 
1 1 2,g 



C 



L,g I 



'19 



c 



JV(1) 



_ (7 JV(1) w( ) , ^L,g C N{1) AT(O) , 2/ 3 



c 



AT(2) 



WV(1) 



rfiv(i) 



(45) 



It is evident that the representation in terms of Mellin-moments is of advantage when compared 
to the corresponding x-space expressions. In the latter case one has to find the inverse Mellin 
transforms of quantities where Cf k and 7^ appear in the denominators of the expressions above 
which in general is not possible. 



3.2 F 2 (x, Q 2 ) and dF 2 (x, Q 2 )/dt 

A second example concerns the structure function F 2 (x,Q 2 ) and its slope. Both quantities 
are well measurable in the present-day deep inelastic scattering experiments. The observables 
Fa,b ( x , Q 2 ) are here 



f b(Q 2 ) = p2 (s)N (Q 2 ) 



This example has been considered before in p7| . In leading order one obtains 

K, 

1 



-^22 



K. 



N(0) 
d2 



JV(0) JV(0) _ N(Q) JV(0) 

iqq '99 iqg 'gq 



N(0) 
2d 



K 



AT(O) 
dd 



-4 

N{0) 



7 



qq 



~ '99 



(46) 



(47) 



The next-to-leading order kernels read : 

: 



K, 

N(l) 



22 

JV(1) 
l 2d 



K. 



d.2 





1 

4 



N(0) N(l) + 

'99 1 qq j "" 



gg 'qq 



N(l) N(0) 

iqg 'gq 



N(0) JV(1) 

iqg 'gq 



(48) 
(49) 



J^L ^(0)^(0) _ ^.N(0)M0)\ , Po r N(l) ( JV(0) 



2/3 



iqq 'gg 'm 'qg 1 n 2 >q I 'aa 



• <)<) 



1 gq 



> qg 



N(0) 

gg 



2/3 
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K 



JV(1) 
dd 



n r N(l) 

" 2 7g?° } 



( AW 



'99 



^(O)^JV(O) , 9/v JV(0) JV(0) 
~ '99 '99 



7 g q 7. 



.9.9 



2/3o 7( Ar{0) 



V JV(0)^(1)' 



A) ( AT(l) _ Tgg %g 
o I '99 AT(0) 
J V 799 



(50) 



N(l) , JV(1) _ A / JV(0) , V( ) 
9 A) V 99 



79 V 9 (0) 



^JV(l) 



V '99 '; 



2V(0) 
99 



2A)-7£ W 



+4/5 c 2 a ; (i) -2a 



(51) 



For this combination in next-to-leading order the evolution depends on two evolution kernels 
only. In the case of polarized deep inelastic scattering similar relations apply considering the 
structure function g 1 (x,Q 2 ) and its slope. The anomalous dimensions and coefficient functions 
of the unpolarized case have to be substituted by those for polarized scattering. 



Although enforced by Eq. (|SE| ) one has still to show that the kernels ( 0-5SL f^THST]) are 
scheme-independent by an explicit calculation, which we have done using Eqs. (^,0) for the 
next-to-leading order contributions. In leading order the scheme-invariance is visible explicitly, 
since the leading order anomalous dimensions and the lowest order coefficient functions for 
Fl{x,Q 2 ) are scheme invariants. 



With the help of the evolution equation (|3q) we are now prepared to ask for the validity of 
crossing relations between different space- and time-like quantities in perturbative QCD. Rela- 
tions of this kind are henceforth called Drell-Levy-Yan (DLY) relations, although the original 
reasoning of these and other authors was quite different. One condition to ask such a question 
at all is that the behavior of all contributing parts under crossing from space- to time-like 
momentum transfer are controlled. At large momentum transfers \q 2 \ the single parton picture 
applies and the non-perturbative parton densities factorize. This makes it possible to study the 
respective evolutions kernels without reference to the non-perturbative input densities. Even if 
a crossing relation for these quantities does not exist, one still may investigate whether it exists 
for the perturbative evolution kernels. A further condition for this investigation is that the latter 
quantities are scheme-invariant, as in Eq. (|36l). 



4 Drell-Levy-Yan relations 

In the following we study in detail an interesting relation between deep inelastic lepton hadron 
scattering and e + e~ annihilation into a hadron and anything else, proposed by Drell, Levy and 
Yan 0. Here we first briefly illustrate the idea behind the work of DLY for completeness. In 
field theory, the deep inelastic e + e~ annihilation can be related to matrix elements of hadronic 
electromagnetic current operators similar to that of deep inelastic lepton-hadron scattering. 
The crucial difference, apart from the ones which originate from the kinematics, is that the 
annihilation process is not related to the forward Compton amplitude contrary to deep inelastic 
scattering because in the former process the hadron is observed in the final state. Nevertheless, 
both processes are related by crossing symmetry which any field theory enjoys. This motivated 
DLY to study the process in detail and then relate it to the deep inelastic scattering process. 
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From the structure of the hadronic tensors W^ u (q,p) (space-like) and W^(q,p) (time-like) and 
using the standard reduction formalism one can infer that 

W^( q ,p) = -W* u (q,-P) , (52) 

where the momenta within the respective parentheses of the above quantities are the same as 
those defined in the beginning of the paper. In the Bjorken limit for both deep inelastic scattering 
and deep inelastic annihilation for q 2 = —Q 2 ,p.q — > oo and q 2 = Q 2 ,p.q — > oo, respectively the 
scaling structure and fragmentation functions satisfy the following relation []: 

F l {S \x B ) = -(-l) 2 ^+^x E FP(—) , f = l,2,L. (53) 

Here it has been assumed that non-perturbative input parton densities can be decoupled trivially 
and are the same. In other words, the functions F^ t \xe) are the analytic continuations of the 
corresponding functions F> (xg) from 0<£B<ltol<££<oo. This is true only when the 
continuation is smooth, i.e. if there are no singularities for example at x — 1 etc. This relation 
is called DLY-relation in the literature. 

In this section, we study this property in more detail extending earlier work JET]. It is par- 
ticularly interesting to study the above transformation at the level of the splitting functions and 
coefficient functions which constitute the physical quantities such as the structure and fragmen- 
tation functions. Then we show how these relations are preserved for the physical quantities by 
looking at the kernels discussed in the previous section. Apart from scaling violation one also 
encounters distributions of the type 

8(l-z) , ( ln i (1 _~ Z) ) , 1 = 0,1,2,..., (54) 

which destroy the continuation through z = 1. Here the distribution (ln*(l — z)/(l — z)) + is 
represented by 

\ 1-z J + (z + l) 0--Z) 

where 5 <C 1. It turns out that the DLY-relation is violated for the coefficient functions and 
splitting functions separately because both are scheme dependent. This in particular happens 
when we adopt the MS-scheme. Here the relation is already violated up to one- loop order for 
the coefficient functions. Although one can choose other schemes in which Eq. ( |5*3"| ) is preserved 
P"7jj ) up to one-loop order we do not know whether this will hold up to any arbitrary order 



sec 



in perturbation theory. 

Let us start with the simplest examples and consider the scheme-invariant evolution kernels 

Eq. 



4.1 The Drell-Levy-Yan Relations at Leading Order 

In the case of the scheme-independent evolution kernels describing the evolution of F2(x, Q 2 ) 
and dF 2 (x, Q 2 )/dt, respectively, or its polarized counterpart for the structure function gi(x, Q 2 ), 

5 Here we indicate the overall signs in case of the scattering of particles of different spin, cf. |30|. In the 
original work of DLY M the Yukawa-theory was discussed which does not contain gauge bosons. 
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only the transformation of two combinations of the leading order anomalous dimensions has to be 
considered, cf. (f|7|)- These are the determinant and the trace of the singlet anomalous dimension 
matrix at leading order. In both quantities the color factors of the off-diagonal elements enter 
only as a product. The unpolarized and polarized leading order splitting functions read 



p(0) 

qg 



AP (o) 
qg 



p(0) 

gq 



AP (o) 

gq 



p(0) 

gg 



AP (o) 
gg 



4C7 



8T R N f z 2 + (l- z f 



z 2 -(l-z? 



8T R N f ^ 

l + (l-z) 2 



4CV 



8CL 



l-(l-z) 



z 
z 



8C 



A 



l-z) 



l-z) 



1 - z ft \ 
+ + zil - z) 



+ 2P 6(1 - z) 



1 - 2z 



2Po6{l-z) 



(56) 

(57) 
(58) 

(59) 
(60) 
(61) 

(62) 



The crossing relations of the leading order splitting functions are 



p(0) 

qq 



p(0) 

gq 



-zP® - 



qq 



MfT f 

C p 



p(0) f I 
gq 1 



p(0) 
9S 



p(0) 
ff9 



c F 

2 .X I T qi > 



(63) 



where one demands 



6(1 -z) ->-<J(l-z) 



(64) 



Eq. (|63D is easily verified and implies the validity of the crossing relation from space- to time- 
like evolution kernels Eq. (^), i.e. the validity of the DLY-relation for this case. For the 
second set of physical evolution kernels the DLY-relation follows at leading order referring to 



the transformation relations for the leading order longitudinal coefficient functions, Eqs. (|7q , [79]) 
in an analogous way. 



4.2 NLO Splitting function 

As we know, the splitting functions and coefficient functions are not physical quantities due to 
their factorization scheme dependence. Hence, the naive continuation rule for these quantities 
may be violated, which is indeed the case in most of the schemes, e.g. in the MS scheme 
characteristic of n-dimensional regularization. It was demonstrated by Curci, Furmanski and 



Petronzio |33] that by an appropriate modification of the continuation rule in the MS scheme 
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one can show that the time-like splitting functions are related to their space-like counter parts. 
Since the modification of the continuation rule has to do with the scheme one adopts, it simply 
amounts to finding finite renormalization factors. It was shown that the finite renormalization 
factors can be constructed from the e-dependent part of the splitting function when computed 
in dimensional regularization [34]]. In addition to this, care should be taken when dealing with 
quark and gluon states which was not the case in the work by DLY, where a color and flavor 
neutral field theory was discussed. The transformation rules are : 

• The diagonal elements of the space-like flavor singlet splitting functions P qq , P gg have to 
be multiplied by (—1). 

• The off-diagonal elements of the singlet splitting functions matrix have to be multiplied 
by Cp I \2NfTf) for P qg and 2NfTf/Cp for P gq , respectively, accounting for the interchange 
of the initial and final state particles under crossing. 

Note that these transformations are automatically accounted for in the case of the leading order 
physical evolution kernels discussed in the previous paragraph. Keeping this in mind and using 
the known splitting functions OTL |33], p5l and the continuation rules 





ln(l - 


z) -> ln(l 


— z) — ln(z) + in , 




(65) 




\n(S) - 


-> ln(5) + in , 




(66) 


one finds that 












p(l)(S) _ p(l)CT) _ 
qq qq 


-2/?o^ )(1) 








(67) 


p(X)S _ p(l)(T) _ 
qg gq 


-2/5o^ )(1) 


+ zj?w s 


(IS? - *$>) + PS? « 


' (z£ )m 


-ZgW)) j( 68) 


p(l)S _ p(l)(T) _ 

gq qg 


-2/3o^ )(1) 




<3? - 3?) + 3? « 




-4PW),(69) 


pO)S _ p(l)(T) _ 

gg gg 


-2/?o^ )(1) 


+ *™ « 






(70) 


where the quantities with 


a bar denote 


that they are continued from z — > 


1/z with 


the appropriate 



factors in front. These quantities read in explicit form 



^(z) 



P 9q (z) (n) 



. 7 pW 
^ r qq 



2N f T f 



4^) 



(n) 



Cj? z p(n) 

2iV / T / 99 



■4 n) ' - 



(71) 



The relations given in Eqs. (|67|-7T) remain true for the polarized splitting functions []32], [34] as 
well. The renormalization factors appearing in the Eqs. (|6"T|-|7Tj|) are given by 



7 T(1) 



d(0) 



( ln(z) + dji 



(72) 



The constants are different in the unpolarized and polarized case. For unpolarized scattering 
they read 

, d qg = — , CL gq = — , (Jo) 



"qq 



' gg 



1 

2 



'gq 



1 

2 ' 
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whereas in the polarized case 

(Hi = • (74) 

The logarithms in the renormalization factors originate from the kinematics. In dimensional 
regularization, when one continues the partonic structure function J 7 ^ in Eq. (|l^) from the 
space-like to the time-like region one obtains an additional factor z e which when multiplied 
with the pole in e yields ln(z). Since the pole is always associated with the splitting functions, 
one has the function P$' along with ln(z). The ^-independent constant a^, which is also 
multiplied by the splitting function, results from the polarization average. For deep inelastic 
scattering one averages the processes with one gluon in the initial state by a factor l/(e + 2). 
Such an average is not needed for the annihilation process since here the gluon appears in the 
final state. Notice that the average over the polarization sum does not show up in the polarized 
structure functions. Therefore in this case the constants dy are zero. 

The transformation behavior of the non-singlet splitting functions in NLO have been worked 
out in where also the relations for the NLO non-singlet coefficient functions were presented. 



4.3 NLO Coefficient Functions 



Now, let us study how space-like and time-like coefficient functions are related. The coefficient 
functions are expected to violate the DLY-relation due to their scheme dependence. Here we 
first present the relations between the space-like and time-like coefficient functions C^{z) (i = 
l,L;k = q, g). The leading order transverse coefficient functions are identical. At next-to-leading 
order, in the MS scheme ||36|| , the coefficient functions are related by the Z-factors in Eq. ( [72] ) 
as follows : 



(T)(l) 
,9 1 



c 



2N f T, 



{z)+[zC[ 
2zC{ 



(S)(l) 



f 1 ! 



,9 



(75) 
(76) 



Since the coefficient functions depend on the hard scale of the process, one has to replace the 
space-like q 2 by the time-like q 2 in addition to Eqs. (|65|- |66|) . This leads to the following contin- 
uation rule 




space— like 




— ITT 



(77) 



time— like 



The Z-factors get contributions from two sources. The first one is ^-dependent and comes from 
the phase space integrals. The time-like phase space acquires an extra factor z e which gives a 
finite contribution when being multiplied with the pole terms 1/e. The pole term originates from 
the collinear divergence in n-dimensional regularization. The second term originates from the 
polarization average which is again absent in the time-like case. The continuation rules given in 
Eqs. (|65]-|6^, [77] ) are essential to get the constant £(2) right when one goes from the space-like to 
the time-like region. Notice that the space-like coefficient function contains —4^(2)5(1 — z) and 
the time-like one contains 8£(2)5(1 — z). The difference which is 12£(2) can be understood to 
originate from the one-loop vertex correction when one continues from the space-like to time-like 
region in Q 2 . The same also holds when other regularization methods for the collinear divergences 
are chosen. It is worth noticing that if one would replace ln(l — z) — > ln(l — z) — ln(z) contrary 
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]) one would obtain an additional term 12£(2) on the righthand side 



to the prescription in Eq. 
of Eq. fl7§). 

The zeroth order longitudinal coefficient functions are identically zero so that the first order 
contributions are scheme independent. This implies that there are no pole terms in the corre- 
sponding partonic structure function Tl^- Hence, there is no left-over finite piece which could 
arise from the z e - or n— dimensional polarization average. We find 



C 



S (1) w 

C F 



2 







2N f T t 

4.4 NNLO Coefficient Functions 



zC 



(5)(1) 
L,g 



(78) 
(79) 



4.4.1 Longitudinal Coefficient Functions 



We consider the NNLO correction to the longitudinal coefficient function. We follow the results 
given in [B7|, |D| for the space-like and [fFJ, for the time-like case. It turns out that the 
coefficient functions are related by the Z-factors through the matrix-valued convolutions 



-C 



(S)(2) 
L,q 
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2N f T f 



z r (s)(i) 
2 ,a 



30) 



< K2) (^) + 
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2N f T, 



f 1 ! 
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(S) (2) 
L,g 



<?9 



Z r (S)(l) 

2 



+Z (T)(1) 
^ 99 



2NfT, 



f 1 / 



-C 



(1)5 
L,9 



ID 



The right hand side of the above equation contains the convolutions of Z-factors with the 
continued NLO longitudinal space-like coefficient functions. We have found this pattern by 
comparing the scheme transformation which we derived in the last section. The reason for this 
structure relies on the fact that Cx $ is obtained as the difference between C 2) ; and Ci j. Since 
the NLO coefficient functions involve various Nielsen-integrals, we used the following identities 
to simplify the expressions : 



Li 2 



Li 2 (-z)-lln 2 (^)-C(2) 



(82) 



Li 2 1 - - 

z 



iln 2 (^)-Li 2 (l 



(83) 



Si, 



1-i 

z 



iin 3 (z) + S 1)2 (l-z), 
o 



34) 
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Li, 1 



6 



\n 3 {z) + Si, 2 (l -z)- Li 3 (l - z) + ln(z)Li 2 (l - z) , 



Lis - 



s,(-i) 



Li 3 (-z) + -ln 3 (z)+C(2) ln(z) 
o 



-S 1>2 (-z) + Li 3 (-z) - ln(z)Li 2 (-z) - - ln 3 (z) + C(3) 

o 



i6) 



(87) 



If we do not continue ln(l — z) and replace ln(l — z) — > ln(l — z) — ln(z), then terms propor- 
tional to C(2) are not compensated between space-like and time-like coefficient functions and 
hence the relations given in Eqs. (j50| , 31) are no longer true. Although formally of NNLO, the 
coefficient functions C^U (z) may be combined to physical evolution kernels together with the 
NLO splitting functions as shown in section 3.1. In Section 4.5. we will show that because of 
the transformation in Eqs. fl8"0| , |B~TD the physical evolution kernels in sections 3.1 and 3.2 remain 
DLY-invariant . 



4.4.2 Transverse Coefficient Functions 



In NNLO physical evolution kernels for the transverse structure and fragmentation function can 
only be constructed when the space-like and time-like three-loop splitting functions are known. 
If they become available one can extend Eqs. ((47H5T|) up to second order. Here we consider the 
relation between the space- and time-like coefficient functions using the transformation relations 
(|6"4] , |65| , p6|) and ([T7|) for unpolarized and polarized scattering. 

The space-like coefficient functions for unpolarized scattering are computed in []39], |40 | 
whereas the time-like ones can be found in [13|]. The transverse coefficient functions are 
related by (see appendix 6.1) : 



C^\z) + {zC[T ) (± 
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~ qq 

+ Z (T)(1) 

' gq 



p(0) 

gq 



Hz) + \z™ zg>n 



C F r (S){i) 



+12C|C(2) 2 In 



qg 



-3 5(1 -z) 



For the polarized NNLO coefficient functions which were derived in [ p9| , |4"0|j and |42|, [43|, we find 
that the form of Eqs. ([88]) is the same but the term 



does not occur. 

Similarly for the gluonic coefficient functions we find 
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(S)(l) 



5*^ + ^®^"^) ' < 89 > 



For the polarized case, the terms 



1/3 p(o) + J_p(o) 

99 16 99 



3?) 



in Eq. (p9|) are absent. Since we do not have to average over the initial state gluon polarization 
in the case of polarized scattering the term hx(z) + 1/2 multiplying the first bracket in Eq. (| 
is replaced by ln(z), cf. also Eq. ([F3, [T4]). 



4.5 NLO Physical Evolution Kernels 

After having found the relations between space-like and time-like splitting and coefficient func- 
tions, we investigate the DLY-transformation for the physical evolution kernels presented in 
sections 3.1 and 3.2. In order to do this, we define the difference between the time-like quanti- 
ties Kfj and the continued space-like quantities K^ by 
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(90) 



where Kfj is obtained by transforming Kfj to the time-like region using the continuation 
rules (65], |66| , |64| , 77) . Application of the DLY-transformations provides us with the follow- 
ing results 



SK. 



JV(1) 
22 



^L,g r jy(i) ^L,g c C N(l) N(0 ) 
^N(1) U Igq ~ ^yNO) 2 >9 Iqq 



^L,g 



+ 



r N(l) ■ 
r N(l) 



sr N(l) 



^iV(l) r r N(2) 

^JV(l) ^AT(l) 
°i,9 U i,9 



99 



c 



JV(1) 



c 



7V(1) 



U L,g 



C 



NO) 2 >9 
i,9 



_ oa v(T)iV(l) _ ~JV(0)^(T)JV(1) , ^AT(O) v(T)iV(l) 
/q<j rO^oo '00 00 '00 00 



99 



'99 



' 99 



.99 



c 



AT(1) 
L,o 



r»;v(i) 



_L na 7 (T)N(l) _ 7 (T)7V(1) -7V(0) , 7 (T)7V(1) -7V(0) 
'99 ^r^O-^qq qq Iqq ' qq Iqq 



9.9 



9.9 



99 



00 



' 9,9 



_ v(T)7V(l)-/V(0) , Z (T)NO)^N(0) 
gg "oo ^ oo 'oo 



' 9.9 



9.9 



' .9.9 



(91) 



19 



Substituting the expressions for and Sj^ 1 ^ using Eqs. (|67| , |69|) , we get 

SKgV = . (92) 
For the remaining evolution kernels one obtains 
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The explicit expressions for the differences in the coefficient function are given in appendix 6.1 as 
well as a series of involved Mellin-convolutions leading to Nielsen-integrals (see appendix 6.2), 
which are necessary in the explicit calculation. 



Using Eqs. (p7| -|70f), leads to 
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The physical evolution kernels Kj j for the evolution of the structure functions F 2 and Fi are 
thus DLY-invariant to next-to-leading order if continued from the space-like to the time-like 
region. 

We turn now to the physical evolution kernels in next-to-leading order where we choose the 
physical quantities F 2 , dF2/dt as a basis. Here only two evolution kernels are contributing, which 
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change under the DLY-transformation as follows : 



SK. 



iJ2 



| (< W " Z^A + t2<°> - 2/3 



2 -N(0) y 0L, 2g ^qg J \[l qq ) 1 qq 1 gg 



lqq 



(99) 



SK, 



0o 



dd 



Z (T)N(1) 
99 



, '99 



^(0) 
99 



2A 



+4/? (<«-^«). 
From Eqs. ( |75| , |7B| , |75| , |79|) we can derive that 

5^2 = 0, 



SK 



dd 



. 



(100) 



(101) 
(102) 



From these results it is clear that the time-like physical evolution kernels K^ can be directly 



derived from the space-like physical evolution kernels using the continuations in Eqs. (B 



51 , |77| ) where one has to account for the corresponding changes in the overall color factors. The 
Z T -factors which are needed for the transformation of the splitting and coefficient functions 
cancel in the expression above. In the future one can extend the investigation performed in 
this section to physical evolution kernels at the NNLO-level, provided the 3-loop anomalous 
dimensions are calculated. For the choice of observables (F 2 ,F L ) one also needs the three-loop 
coefficient functions. 



We finally would like to comment on a relation derived by Gribov and Lipatov in [19| for the 



leading order kernels for a pseudoscalar and a vector field theory. [] One may write it in the form 

K(x E ,Q 2 ) = K(x B ,Q 2 ) , (103) 

where K and K denote the time- and space-like evolution kernels, respectively, and xb = 1/^e- 
One verifies, that this relation holds in leading order for the space- and time-like splitting 
functions of QCD, Eqs. (|56|-[62|), without changing the ^-function, Eq. (|64|). 

Starting with next-to-leading order, this relation is not preserved. For the physical non- 
singlet evolution kernels this was shown in ||33 , j35fl and for some singlet combinations in [|35|| . We 
find, that also for the physical singlet combinations, Eqs. (f|2|-f|5], |48[-|5TD, this relation is violated 
as well. 



5 Conclusions 



The old question, whether the scattering cross sections of deep inelastic scattering e +P — > e + 
l X' are related to the annihilation cross section e + + e~ — > P+ l X' by a crossing relation changing 

6 See also fEil for related work. 
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from t- to s-channel was newly discussed. Since in both reactions non-perturbative quantities 
such as the structure and fragmentation functions contribute the above question cannot be 
answered by means of perturbation theory for the process as a whole. However, since both the 
parton densities involved in the space- and time-like process factorize if the virtuality Q 2 = 
\q 2 \ of the four-momentum transfer is large a related question can be asked for the crossing 
behavior of the respective evolution kernels, which are computable within perturbation theory. 
In the calculation of both inclusive processes only two types of singularities occur, the collinear 
singularity and the ultraviolet singularity. These divergences are absorbed into the bare parton 
densities and the coupling constant, respectively. Two distinct renormalization group equations 
are implied. They quantify the impact of the factorization and the renormalization scale on the 
DIS structure functions and fragmentation functions when the perturbation series is truncated 
up to a given order. However one can construct factorization-scale independent evolution kernels 
which describe the scheme-invariant evolution of these physical quantities in terms of a kinematic 
variable given by Q 2 . This scheme invariant evolution is guaranteed up to any finite order in 
perturbation theory. Notice that in finite order this method does not remove the dependence 
of the physical quantities on the renormalization scheme of the strong coupling constant or its 
scale \i r . 

The first example of the application of the physical evolution kernels is the coupled structure 
functions F 2 (x, Q 2 ) and Fl(x, Q 2 ) associated with the corresponding fragmentation functions in 
e + e~-annihilation. A second example is given by F 2 (x,Q 2 ) and dF 2 (x, Q 2 )/dln(Q 2 ). Contrary 
to the splitting functions (anomalous dimensions) and coefficient functions the evolution kernels 
of the examples above are factorization scheme independent. For that purpose transformation 
relations have been derived for the splitting functions up to NLO and the coefficient functions 
up to NNLO. We have also shown that these kernels are invariant under the Drell-Levy-Yan- 
transformation up to next-to-leading order. On the other hand the Gribov-Lipatov relation, 
which is valid in leading order, is already violated at next-to-leading order. It remains to be 
seen how the physical evolution kernels behave under the DLY crossing relation at NNLO, which 
presumes the knowledge of the yet unknown three-loop splitting functions (space- and time-like) 
as well the three-loop longitudinal coefficient functions in the first example above. 
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6 Appendix 

6.1 Coefficient Functions 

In this appendix we list the difference of the space- and time-like coefficient functions in the 
MS scheme, which are used in section |] to study the validity of the DLY-relation. Here the 
expressions also contain the logarithms 

L N = ln(Q 2 //iJ) (104) 

which arise when the factorization scale /i| is chosen to be different from Q 2 . 

The difference between the longitudinal non-singlet coefficient functions corresponding to 
the processes l* + q^q + g + g and 7* — > 'q' + q + g + g respectively are given by 



A{2z - \n(z)) ln(z) - 16Li 2 (l -z) + 8-8z] (105) 
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where L q' denotes the quark in the final state which undergoes fragmentation into a hadron 
P (see Eq. (8)). The difference between the longitudinal purely singlet coefficient functions 
corresponding to the processes J* + q — > q + q + q and 7* — > £ <f + q + q + q respectively are given 
by 

SCj^ ps = N f T f C F -8 ^6 + 4z - ^z 2 ) ln(z) + 16 ln 2 (z) - 16 



304 



9z 



+ 64z 



128 
IT' 



(106) 



The same is done for the longitudinal gluonic coefficient functions corresponding to the processes 
7* + g — > g + q + q and 7* — > t g' + g + q + q, respectively. The difference in the coefficient function 
is given by 



2 \ 24 
1 + - - z ln(z) + 8 ln 2 (z) - 28 + — + 



z J 



4z 



-CaCf 



I (j ( 4 - - + z - ^z 2 ^j ln(z) - 16 + ln 2 (z) 



-32(l--)Li 2 (l-z) 



248 
~9z~ 



40 , 

- 24z + — z z 



9 



(107) 



Notice that for the computation of the coefficients functions above and the ones following here- 
after one also needs the virtual contributions to the zeroth and first order partonic processes. 
The differences between the transverse coefficient functions emerge from the same processes as 
mentioned above Eqs. ( |105| , |106| , 107|) . In the non-singlet case we have 



r r (2)NS 
0U l,q 



1 

— I 

2 



Hz) 

'l + z 



2C(2) - 41n(z) ln(l + z) + ln 2 (z) 



+ 4 2(1- z)C(2) + 4(1 - z) + 4(1 - z) ln(z) ln(l + z) 



+2(1 + z) ln(z) - (1 - z) ln 2 (z) ln(z) - 16Li 2 (-z) ln(z) 



x 



l + z 



-l + z 



+NfTfCp 
+CaCf 



8 / 10 



9 V l-z 
M*0 /67 , , 2 



1 + llz ln(z 



4^^ + h» W -2C(2)) + 2 f-f, 



+2(1 + z) ln(z) - (1 + z) ln 2 (z) ln(z) + 4((2)(1 + z) ln(z) 



+C| 



ln(z) 
T-z 



L H . -6 + 4 ln(l - z) ) ln(l - z) + 6^ - 1 
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+ ( - 4L M/ + 6 - y \n(z)) \n(z)j + 4^ - 4(1 + z)L„ f + 1 + bz 

-2(1 + z) ln(l - z)^ ln(z) ln(l - z) + 2 ^ - 2(5 + z)L H + 14 + 40,2 
+2(1 + z) \n(z) ln(l - 2) + 6(1 + z)L^ f \n(z) - 8(2 + z) \n(z) 
+7(1 + 2) ln 2 (^ \n{z) + 2Li 2 (l - z) ^4 + 7 + 7^ ln(z) 



\ / 24 
-10 + 22,2 J +8 (-- — - 



+ 13 + Viz) Si, 2 (l - 2) + 36(1 - 2) 



+32C(2) 



1-2 



1 - z ) ln(z) + 12C(2)(9 - \2L H + 4^)5(1 - z) 



. (108) 



For the purely singlet difference we obtain 

5C[ 2) q PS = iV / T / C F [(-8(4 + 6^+^ 2 )(L M/ +ln(l 

x ln(z) - 160 - — - 112,2 - — z 2 + 4(1 + z) U ln(l - z) + 4L 
9z 9 V 



:)) - 1(3(1- — + 2,2 



10 



38 



+— \n(z) \n(z) \n(z) — 8(1 + - z — —z + ln(l - z) 



9z 



9 



+ 16 \2(1 + z) ln(z) - 2 - 3,2 - -z 2 j Li 2 (l - 2) + 32(1 + z)Si )2 (l - z) 



1168 224 640 1808 



9 27,2 9 27 
The difference between the gluonic coefficient functions equals 

*Cg> = 2C A C F [ (l88 + ™ + 662 + ±fz 2 + (48 - y + 16* + f L w 

+ (~ 4 ~ ~3~r + + ( 40 ~ t + 122 + y^ 2 ) ~ ^ 

+4 (2 + 2 - + 2) ln(l + z) - 16 (l + + ,2) ln(z)L M/ 

+8 ^2 + 2 + ^ 111(2) ln(l + z) - 16 (l + ^ + \n(z) ln(l - 2) 

- (f + £ + 1 z ) lAz) + 4 ( 2 - ! - *) ln2(i - 2) ) lnW 

) ln(l - z) + 2 (2 - 2 - 2) ln 2 (l - 2) 



(109) 
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24 



284 



104 



9z 



9 



+ ( -24 + — + Az - —z 2 ) + ( 32 ( 1 - - - 2z ) \n(z 



z 



+16 (2---z) (ln(l -z) + + 16 + ® + 16z + yz 2 ) Li 2 (l - z) 
+4 (2 + - + z^j (l + 2 \n{z))Li 2 (-z) + 16 (-2 + - + z) Li 3 (l - z) 
8- — - lOzJ Si, 2 (l-^) 
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-12 (2 - --z\ ln 2 (l - z)^ ln(z) + (^64 - — - 18z 
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z^j Mn(l - 2) 



6 



6 32 Li 2 (l - z) + 16 2 2 Li 3 (l - 2 



S ( 10 - ^ - 5z\ Si, 2 (l - 2) - 169 + ^ + 50z 



(110) 



We have computed the same differences between the coefficient functions corresponding to the 
structure function g\(x,Q 2 ) which describes polarized scattering. The analogues of Eqs. (|108 , 
109|, |110j) are given by 
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'l + z 



2C(2) - 4 ln(z) ln(l + z) + ln 2 (z) 
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SAC[ 2) q PS = N f T f C F 
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-8(4 + z) ln(z) ln(l - z) - -(26 + 5z) ln 2 (z) - 4(2 - z) 
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6.2 Convolutions 

Here we list the convolutions of a series of functions, which are needed for the investigation of 
the DLY-relation in section f|. Using the definition in Eq. fllBf ) we obtain 
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